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ABSTRACT 

We have derived ages, metallicities and enhanced-element ratios [a/Fe] for a sample 
of 83 early-type galaxies essentially in groups, the field or isolated objects. The stellar 
population properties derived for each galaxy corresponds to the nuclear re/8 aperture 
extraction. The median age found for Es is 5.8 ± 0.6 Gyr and the average metallicity 
is +0.37 ± 0.03 dex. For SOs, the median age is 3.0 ± 0.6 Gyr and [Z/H] = 0.53 ± 
0.04 dex. We compare the distribution of our galaxies in the H/3-[MgFe] diagram with 
Fornax galaxies. Our elliptical galaxies are 3-4 Gyr younger than Es in the Fornax 
cluster. We find that the galaxies he in a plane defined by [Z/H] = 0.99 log o-q - 0.46 
log Age - 1.60, or in linear terms Z oc ctq ■ Age~°'^. More massive (larger (Tq) a-nd older 
galaxies present, on average, large [a/Fe] values, and therefore, must have undergone 
shorter star-formation timescales. Comparing group against field/isolated galaxies, 
it is not clear that environment plays an important role in determining their stellar 
population history. In particular, our isolated galaxies show ages differing in more than 
8 Gyr. Finally we explore our large spectral coverage to derive log(0/H) metallicity 
from the Ha and [Nil] A6584 A and compare it with model-dependent [Z/H]. We find 
that the 0/H abundances are similar for all galaxies, and we can interpret it as if most 
chemical evolution has already finished in these galaxies. 

Key words: galaxies: stellar content - galaxies: abundances - galaxies: elliptical and 
lenticular - galaxies: nuclei - galaxies: evolution 



1 INTRODUCTION 

A decision over the most appropriate scenario of formation 
and evolution of early-type galaxies is still an open ques- 
tion despite the enormous amount of work already put for- 
ward to model primary observational constraints and scal- 
ing relations. Most recent flavours of models ranging from 
monolithic collapse scenarios to hierarchical clustering pre- 
dictions include Arimoto & Yoshii (1987, 1989); Kauffmann, 
White & Guiderdoni 1993; Bruzual & Chariot (1993); Bres- 
san, Chiosi & Fagotto (1994); Worthey (1994); Einsel et 
al. (1995); Navarro et al. (1996); Steinmetz (1996); Tan- 
talo et al. (1996); Bressan, Chiosi & Tantalo (1996); Gib- 
son (1996a,b, 1998); Gibson & Matteucci (1997); Tantalo, 
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Chiosi & Bressan (1998); Chiosi & Carraro (2002); Terlevich 
& Forbes (2002); and references therein. 

Briefly, in the conventional monolithic collapse scenario 
elliptical galaxies formed as a result of a single violent burst 
of star formation at high redshift (possibly z > 3) and 
evolved quiescently ever since. By contrast, in the Hierarchi- 
cal model, there is no intrinsic difference between star for- 
mation occurring at different epochs. AH galaxies are viewed 
as similar star-forming systems in which there is an equilib- 
rium between the inflow of gas and the rate at which it is 
either consumed or driven out of the galaxies by a super- 
nova wind. Elliptical galaxies are formed by merger and/or 
accretion of smaller units over a time scale comparable to 
the Hubble time. 

As pointed out by Trager et al. (2000b, hereafter TOOb), 
the evidence for intermediate-age stellar population (be- 
tween 1 and 10 Gyr) would favour hierarchical models, which 
more naturally have extended star formation over time. On 
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the other hand, WiUis ct al. (2002) have presented an anal- 
ysis of the mean star formation history and space density 
of a sample of luminous field (in their definition, field is 
equivalent to non-cluster) early-type galaxies selected over 
the rcdshift interval 0.3 z ^ 0.6. In all, the mean star 
formation history of the sample is characterised by an ap- 
parently old (zy > 1), solar to slightly above-solax metaJlic- 
ity luminosity-weighted stellar population that has evolved 
passively since the formation epoch. The mean properties of 
Willis ct al. sample arc markedly similar to the properties of 
morphologically-selected luminous elliptical galaxies in rich 
cluster environments at redshifts z < 1 (Ferreras et al. 1999). 
Though neither result in isolation constrains the extent to 
which early-type galaxies in field or cluster environments 
represent a co-eval or co-metal population, the broad sim- 
ilarity between the star formation history of the dominant 
stellar mass component in such galaxies is consistent with 
similar formation conditions for each population. 

Furthermore, the recent results from the K20 survey 
(Cimatti et al. 2002) have identified a class of old passively 
evolving objects (old-EROs) with derived minimum age of 
^ 3 Gyr, corresponding to a formation redshift of Zf ^ 2.4. 
Pure luminosity evolution models with such formation red- 
shifts well reproduce the density of old EROs (consistent 
with being passively evolving ellipticals), whereas the pre- 
dictions of the current hierarchical merging models are lower 
than the observed densities by largo factors (up to an order 
of magnitude). In this survey, a population of dusty star- 
forming galaxies is also identified at 0.7 < z < 1.5, equally 
populated as the old-EROs ensemble, and suggests a signifi- 
cant contribution to the cosmic star-formation density at Zf 
~ 1. 

On the theoretical side, Chiosi & Carraro (2002) nicely 
summarized the current situation by highlighting that the 
hierarchical scheme melting together subunits made of gas 
and stars is not the dominant one by which elliptical galaxies 
are made unless it has taken place in the very remote past. 

Independent of the adopted scenario for the formation 
and evolution of elliptical galaxies, models should meet some 
observational constraints on the history of star-formation of 
these objects. These constraints include ages, metallicities 
and degree of enhancement in a-elements, derived by means 
of line-strength indices and Single Stellar Population (SSP) 
models (e.g., Gonzalez 1993; Worthey 1994; Kuntschner 
2000, 2001; Trager et al. 2000a; TOOb; Thomas, Maraston & 
Bender 2003, hereafter TMB03; Thomas, Maraston & Korn 
2004). 

This paper explores the central stellar populations of 
a sample of local early-type galaxies and search for corre- 
lations among them and with structural parameters. Many 
previous works have studied such correlations (e.g. Tantalo, 
Chiosi & Bressan 1998; J0rgensen 1999; Trager et al. 2000a; 
TOOb; Kuntschner 2000; Terlevich & Forbes 2002). We pay 
special attention to TOOb, as many of the analysis presented 
here were first explored by these authors. We would like to 
highlight the high mean signal-to-noise of our sample and 
the use of the TMB03 set of SSP models for the analysis 
in this paper. TMB03 presented an empirical calibration for 
the synthetic Lick indices of SSP models that for the first 
time extended up to solar metallicity. The new metallicity 
range covered approaches the regime that is relevant for the 
interpretation of the integrated spectra of elliptical galaxies. 



In Section 2 we briefly describe the sample and stel- 
lar population model applied. In Section 3 we explore the 
line-strength Lick indices and the quality of the model cal- 
ibrations. The determination of ages and metallicities, and 
the cnhanced-clement ratio [a/Fe] are presented in Sections 
4 and 5, respectively. We investigate the existence of the 
age-metallicity-cro-[a/Fe] hyperplane for early-type galaxies 
in Section 6. Discussions and conclusions from the nuclear 
stellar population analysis are presented in Sections 7 and 8. 
Finally, in Section 9 wo explore our large spectral coverage 
to compare metallicities derived from different elements, in 
the form of a log(0/H) vs [Z/H] plot. 



2 DATA AND ADOPTED MODEL 

In the present work we will be analysing the spectra for 83 
galaxies with signal-to-noise (S/N) ratio greater than 15 (per 
resolution element FWHM = 6 A) for the re/8 central aper- 
ture extraction. The observations were performed at the 2.12 
m telescope of the Observatorio Astrofisico Guillermo Haro 
(OAGH), at Cananea, Mexico. Eight galaxies in our sam- 
ple were carefully classified as isolated (typical TuUy 1987 
density of 0.08); the sample also hcis 18 galaxies in the field, 
and 57 group members. The group classification is taken 
from Garcia (1993). The sample also contains throe Virgo 
cluster galaxies, namely NGC 4365, NGC 4374, NGC 4754, 
with the respective TuUy (1987) catalogue density of 2.93, 
3.99 and 2.62. In total, the sample consists of 52 elliptical 
galaxies and 31 bulges of SOs or early-type spirals. Note the 
original sample had 86 galaxies, but for the present analysis 
the galaxies NGC 3139, NGC 5854 and NGC 5869 were not 
used because they do not have spectral information at A 
5100 A, and so no ages and metallicities could be derived. 
Therefore, the total sample with age and [Z/H] information 
has 83 galaxies. The full details of the sample, observations, 
data reduction and index measurements are presented in 
Denicolo et al. (2004, hereafter Paper I). 

The spectra cover a wavelength range from ~ 3850 A 
to ~ 6700 A. The median S/N of the sample is 40 per reso- 
lution element. The Lick/IDS line-strength index measure- 
ments were performed separately in each galaxy frame and 
later averaged, the error of the mean was computed from the 
repeated measurements per galaxy (on average, there are 8 
spectral frames per galaxy) . 

Kinematic information for the galaxies was carefully de- 
rived in Paper I. 25 Lick spectral indices were measured and 
calibrated to the Lick/IDS system. Most of the galaxies in 
the sample had the Balmer indices corrected for emission us- 
ing Hq detection. The emission correction with Hq should 
be loss uncertain than the empirical emission correction of 
the Balmer lines using [O in] A5007 A. 

In this work we refer to our sample as the OAGH sam- 
ple, and use the following definitions for the special spectral 
indices: <Fc> = (Fe4383-hFe5270-hFe5335)/3 and [MgFe] = 
[Mgh ■ (0.72 • Fe5270 + 0.28 ■ Fe5335)]l/^ as in Kuntschner 
(1998) and TMB03, respectively. 

TMB03 have constructed stellar population models for 
various and well-defined element abundance ratios. The 
most important ratio is a/Fe, which is the ratio of the so- 
called a-elements (O, Ne, Mg, Si, S, Ar, Ca) to the Fe-peak 
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elements (Cr, Mn, Fe, Co, Ni, Cu, Zn), because it carries in- 
formation on the formation time-scale of stellar populations. 
The whole set of Lick indices of SSP models is presented with 
[a/Fe] ratios of 0.0, 0.3, 0.5. The impact from the element 
abundance changes on the absorption-line indices arc taken 
from Tripicco & Bell (1995, TB95), using an extension of 
the method introduced by TOOb. The models are calibrated 
with the globular cluster data of Puzia et al. (2002). The 
present models now allow for the unambiguous derivation 
of SSP ages, mctallicities, and element abundance ratios, in 
particular [a/Fe] ratios. For this reason we have adopted 
their models throughout this paper. We note that recent re- 
stilts from the same group of authors (Thomas, Maraston & 
Korn 2004) have extended the models of variable abundance 
ratios to the Balmer lines. The key result is that, unlike H/3, 
the as and H7 indices show a marked dependence on [a/Fe] 
ratio. We will discuss this dependence again in the next Sec- 
tion. 



3 EXPLORING NUCLEAR LINE-STRENGTH 
INDICES 

There are considerable shifts in metallicity and age when 
comparing TMB03 models with G. Wortlicy's SSP mod- 
els (W94). Figure 1 presents a comparison between model 
tracks for the same age and [Fe/H] grid steps. Taking the 
position of the galaxy NGC 3379 in the diagram, we have 
interpolated an age of 10.9 ±i| Gyr and [Fe/H] = 0.30 ±m 
dex using TMB03 models, while W94 models result in 11.4 
±iJ Gyr and [Fe/H] = 0.22 ±g:g| dex. The ages and metal- 
licities were interpolated as will be explained in Section 4. 

The metallicity values derived from interpolation of 
the TMB03 model are generally higher than with the W94 
model (in the case of NGC 3379, TMB03 results +0.08 dex 
higher in metallicity than W94). Likewise, the ages are sys- 
tematically younger in TMB03 by comparison with W94. 

We note that the emission correction (see Paper I, Sec- 
tion 4) used in this work is, in general, larger than in pre- 
vious works. The effect of a higher emission correction, i.e., 
a higher H/3 value, is a decrease in age and consequent in- 
crease in metallicity, due to the bending of the model tracks. 
Hence, one should keep in mind that both the use of new 
model tracks and the larger emission correction are responsi- 
ble for differences in age/metallicity values when comparing 
with the literature. 

Note however, that the differences between the age and 
[Fe/H] from W94 and TMB03 models for NGC 3379 are 
not representative of the uncertainties involved in these de- 
terminations. The uncertainty in the age and [Fe/H] values 
for NGC 3379 is of greater magnitude than the difference 
between the models. We will discuss the errors in age and 
metallicity determinations in Section 4.1. 

3.1 Model vs data 

The aim of the initial paragraphs of this section is to decide 
which combination of indices provides a reliable and well- 
understood determination of the stellar population age and 
metallicity. For that, and to assess the consistency of our 
index measurements in relation to the model predictions we 
present index-index plots which are almost degenerate in 
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Figure 1. Comparison between Worthey (1994, W94) and 
Thomas, Maraston & Bender (2003, TMB03) models. Both mod- 
els are shown in age steps of 1, 2, 3, 5, 8, 12, 15 Gyr and metallicity 
steps of -0.33, 0.00, +0.35. The Wortliey model was interpolated 
to the given age and metallicity steps with the use of G. Worthey's 
home page (http://astro.wsu.edu/~worthey/). The point repre- 
sents NGC 3379. 
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Figure 2. Index vs index plots for Balmcr-linc indices in the 
Lick/IDS system. Overplottcd arc the model predictions by 
TMB03. The solid circles are 83 galaxies from our sample. The 
error bars represent the average error of the mean per index, from 
repeated observations. 

age and metallicity. We separate the indices in three groups 

based upon their sensitivity to a elements, Fe-peak elements, 
and the Balmer lines. This group distinction is also useful to 
check how well the various Lick indices trace the correspond- 
ing element abundance. An analysis of the same type was 
carried out by TMB03, and we feel it is important to show 
the consistency of the models with respect to our sample as 
well. In the following figures we plot 83 galEixies observed at 
OAGH. 

We analyse firstly the Balmer-line indices. Figure 2 
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Figure 3. Index vs index plots for indices that correlate with Mg. 
The symbols are the same as in Figure 2. 



shows the relation between H/3, H7a,f and HJa.f- The H/3 
index is only marginally affected by the [a/Fe] ratio, as dis- 
cussed in TB95 and in Thomas et al. (2004), i.e., more 
important for the modelling is taking correct account of 
the Horizontal Branch morphology (Maraston et al. 2003). 
The models used here have the mass-loss on the red giant 
branch phase taken into account, as calibrated in Maraston 
& Thomas (2000). The agreement between models and ob- 
servations is not as good as one would expect, even given the 
observational errors. However, Thomas et al. (2004) discuss 
about the important dependence of the higher order Balmer 
lines H^,^ on [a/Fe]. To correctly compare the Balmer lines 
then wc would need to previously know the overabundance 
ratio for each galajcy and use the corresponding model grid 
for each galaxy [a/Fe] value. This fact certainly diminishes 
the advantage of higher-order Balmer line indices over H/3 as 
age indicators, as we will discuss in Section 4. We then em- 
phasize to the reader that some mismatch between galaxies 
and models of higher-order Balmer lines at constant [a/Fe] 
is expected, because H^^-^ change significantly for different 
[a/Fe] ratios. 

In Figure 3 the a-element sensitive indices Mgi,2, 
TiOi,2, NaD, CN1.2 seem consistent with each other. The 
"Ca lines", Ca4227 and Ca4455, show a mismatch with 
the models; this mismatch grows with increasing Mgb. In- 
deed the model predictions of the Ca4227 index seem to 
be slightly too high as discussed in Maraston et al. (2003); 
those authors indicate that models with variable a/Ca could 
explain the difference. 
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Figure 4. Index vs index plots for Fe-peaJj indices. The symbols 
are the same as in Figure 2. 



We have plotted the Ca4455 index in the same dia- 
gram of the a-element sensitive indices to exemplify how 
misleading can be the actual index nomenclature. Despite 
its name Ca4455 is insensitive to Ca abundance (TB95), 
while Fe and Cr, both elements of the "depressed" group, 
are the dominant contributors to this index. Concerning the 
model tracks, Ca4455 responds strongly to [a/Fe] ratio, and 
it most likely suffers from calibration problems (Maraston 
et al. 2003). 

The TiOg agrees reasonably well with Mg6 at different 
metallicities. NaD also behaves consistently with Mg6, but 
the high sensitivity of the NaD index to interstellar absorp- 
tion severely hampers its usefulness for stellar population 
studies (Worthey et al. 1994). We noticed that probably part 
of the scatter in the TiOi index is due to the broad NaD fea- 
ture, once the blue pseudocontinuum of TiOi is very close in 
wavelength to the NaD profile. The CN indices show a good 
match between the galaxy sample and the models. We ob- 
serve a mismatch of models and galaxy data for the G4300 
index. The calibration of the models with the globular clus- 
ter data for this index is also not convincing as already noted 
by Puzia et al. (2002). 

In the diagrams with the indices most sensitive to Fe 
(Figure 4), we draw attention to the relatively large scat- 
ter in the index values for Fe5709, Fe5782, and specially for 
Fc5015, the last one mainly due to uncertainties in the emis- 
sion contamination correction (the index Fe5015 was emis- 
sion corrected as explained in Paper I, Section 4.2.2). We 
plotted the C24668 index together with the Fe group as in 
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Figure 5. Age-metallicity diagnostic diagrams. The Balmer-line indices are corrected for nebular emission as explained in Paper I. 
Overplotted are the model predictions by TMB03 (with [a/Fe] = 0.0). The dashed lines represent lines of constant metallicity ([Z/H] = 
-1.35, -0.33, 0.00, +0.35, +0.67 dex), whereas the solid lines are lines of constant age (ages = 1, 2, 3, 5, 8, 12 and 15 Gyr). Open symbols 
denote SO galaxies, solid symbols represent elliptical galaxies. Isolated galaxies are plotted with larger symbols. The error bars at the top 
of the plots are the average error of the mean per index. The left-top corner SO galaxy far away from the main distribution of galaxies 
is NGC 3156. 
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Maraston et al. (2003). Although this index is supposed to 
follow carbon, it shows a correlation with Fe. C24668 must 
have calibration problems too and its mismatch with the 
models cannot be explained by overabundance (TMB03). 
As in the globular clusters studies, this index shows a poor 
match between models and data (Puzia et al. 2002). In view 
of modelling uncertainties (TMB03; the mismatch of C24668 
between model and observation cannot be explained only by 
a/Fe effects), C24668 becomes an index not well suited for 
element abundance studies. 

In all, the distribution of the OAGH galaxies in these 
diagrams is very similar to the distribution of the bulk of 
galaxies from the sample of 381 Lick/IDS galaxies in Tragor 
et al. (1998) (see comparison of data and models in TMB03). 

As pointed out by Worthey (1994) the determination of 
the age of old stellar populations is complicated by the sim- 
ilar effects that age and metallicity have on the integrated 
spectral energy distributions. Broad band colours and most 
of the indices are degenerate along the locus of A age ~ - 
3/2 A Z. Worthey dubbed this behaviour the "3/2 rule". In 
other words, the integrated SED of an old (<^ 2 Gyr) stellar 
population looks almost identical when age is doubled and 
Z reduced by a factor of three at the same time. Only a few 
narrow band weak absorption line indices have so far been 
identified which can potentially break this degeneracy. In 
terms of age, the Balmer lines H/3, H7 and US are the most 
promising features, being clearly more sensitive to age than 
to metallicity. Features which are more metal sensitive, and 
thus less age sensitive than the average are: C24668, Mg2, 
Fe5270, Fe5335 and Fe4383. 

In Figure 5 we show metallicity vs age indicators, in the 
form of [MgFe], <Fe>, C24668 against H/3, Hja^f and USA- 
Note that the degeneracy of the models is stronger for older 
metal rich systems, i.e. those in the bottom right corner of 
the diagrams. Both ellipticals and SOs show a spread in age 
of a factor of about 4 (mostly between 3-12 Gyr), and a 
larger spread in metallicitics (approximately from -0.33 to 
-1-0.67 dex). The outlier galaxy on the top- left corner of the 
plots in Figure 5 is NGC 3156, a disk-dominated SO galaxy 
with a rather faint bulge (Michard & Marchal 1994). As we 
have seen, when well calibrated by the models, C24668 could 
be the most promising index for [Fe/H] determinations given 
its high sensitivity to metallicity (its range of variation from 
[Z/H] solar to three times solax is about 5 A). Considering 
ease of measurement, existing calibrations, large metallic- 
ity dependence, small age dependence, and relatively tight 
relation, we can say that, at present, the best metallicity 
indicator in Figure 5 is [MgFe]. Note that we are using the 
[MgFe] index defined in TMB03, where Fe5270 and Fe5335 
are weighted as 72% and 28%, respectively, to balance the 
[a/Fe] dependence (see Fig.7 of TMB03). 

Concerning age indicators, H/3 is still the best calibrated 
Balmer-line index, even though it suffers from emission con- 
tamination. The higher order Balmer lines (i.e., H7, H(5) are 
a lot less sensitive to emission (see Table 6 in Paper I), but 
at the same time they are more affected by abundance ef- 
fects because of the presence of more metallic lines in the 
bluer parts of the spectrum (Thomas et al. 2004). Maras- 
ton et al. (2003) present a calibration of the Balmer lines, 
in U/j^-y^s vs [MgFe] diagrams, for 12 globular clusters with 
[MgFe] ^ 3.5 A. We emphasize that, as in TMB03 (their 
Fig.2, concerning index calibration) the behaviour of glob- 
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Figure 6. The panels compare the ages and metallicities in the 

catalogue by Tcrlcvich & Forbes (2002, TF02) with our data de- 
terminations (OAGH). One-to-one correspondence is indicated by 
the dashed line. We have interpolated ages and metallicities for 15 
galaxies in common with TF02. Both age and metallicity values 
for this comparison were derived from the H/3 vs [MgFe] diagram, 
using Worthey (1994) stellar population models (see Table 1). 



ular clusters does not always represent the range of values 
covered by galaxies. There are significant differences in the 
index calibration for clusters and galaxies specially for in- 
dices such as C24668, CNi,2, Ca4227, Ca4455, NaD, TiOi,2, 
US and H7. The H/3-index calibration for globular clusters 
seems robust, despite the very large scatter of the galgixy 
sample from Trager et al. (1998). 

In this work, ages and metallicities (according to the 
TMB03 models) are derived from the [MgFe] vs H/3 plane 
as shown in Figure 5 (and also in the top panel of Figure 
7). Note that there is some level of mismatch between data 
and models and a significant number of galaxies falls outside 
the models. We interpolated between the model grid points, 
but did not attempt to extrapolate to ages or metallicities 
greater than 15 Gyr or 0.67 dex, respectively (in the plots we 
have denoted these galaxies by arrows and placed them with 
symbolic values for age: 16 Gyr and [Z/H]: 0.69 dex). The 
only galaxy we attempted to extrapolate a consistent age 
and [Z/H] value was NGC 3156. We have used a polynomial 
interpolation of varying order, depending on the position 
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Table 1. Comparison with TF02 of ages and metallicities derived 
using Worthey (1994) models. 



Galaxy Age OAGH Age TF02 [Fe/H] OAGH [Fe/H] TF02 
[Gyr] [Gyr] [dex] [dex] 



ably duo to the emission correction applied in this work. 
Panel b reveals the corresponding small shift to higher metal- 
licities also very likely caused by the higher H/9 values result- 
ing after the emission correction. There is indeed a correla- 
tion between A age (TF02-OAGH) and the applied emission 
correction EC(H/?) (Paper I). 



4 AGES AND METALLICITIES 

The [Z/H] - log Age relation is presented in the bottom 
panel of Figure 7. The small arrows represent lower lim- 
its in age and metallicity for the galaxies that fall outside 
the model grids. Those galaxies for which we can only give 
a lower limit in age are mostly ellipticals, while those for 
which we give a lower limit in metallicity are as many el- 
lipticals as SOs. It is worth emphasizing that the age and 
metallicity values derived are not absolute, only the rank- 
ing of the gala^xies is relatively independent of the adopted 
models. An error ellipse is shown at the bottom corner of 
the plot. This la error ellipse is the result of the interpo- 
lation in the H/3-[MgFe] plane of a galaxy with H/3 error of 
0.11 A and [MgFc] error of 0.06 A, which correspond to the 
mean error values for these indices in the sample. Note that 
the correlation between age and [Z/H] in Figure 7 (roughly 
A log age = —1.57 A [Z/H]) is consistent with being origi- 
nated by correlated errors, as the spread of points has similar 
orientation to the major-axis of the error "ellipse". A Spear- 
man correlation test results in p = -0.566 with more than 
90% probability that a correlation is present. The vectors 
shown in the error ellipse correspond to 0.078 dex in [Z/H] 
and 0.17 dex in log Age; these values are a rough average of 
the overall uncertainties in age and metallicity determina- 
tions for our sample. A further analysis of the errors in the 
age-met allicity diagram is presented in the next Section. 

SOs and Es show a slightly different distribution in the 
plane [Z/H] vs log Age. Excluding the three Virgo clus- 
ter galaxies, the sample of mostly low-density environment 
galaxies is composed of 50 ellipticals and 30 SOs. The median 
age value for Es is 5.8 ± 0.6 Gyr, and for SOs is 3.0 ± 0.6 Gyr, 
where the uncertainties are the errors of the mean. The me- 
dian metallicity value for Es is 0.37 ± 0.03 dex and for SOs, 
0.53 ± 0.04 dex. Therefore SO galaxies seem to be slightly 
younger and more metal-rich than ellipticals. However the 
elliptical galaxies in our sample show a larger spread in age 
distribution than SOs. Note also that the metallicities of the 
young SO galaxies are subject to large errors because they 
are likely to have disk contamination, i.e., composite stellar 
populations. In particular, the isolated galaxies themselves 
also cover a large range in metallicity and age. On the other 
hand, our two Virgo cluster elliptical galaxies are placed rel- 
atively close to each other in age (filled square symbols in 
Figure 7) . It might be expected that cluster galaxies consti- 
tute a more homogeneous sample (as demonstrated in the 
Fornax cluster by Kuntschner 1998), with perhaps a larger 
spread in metallicities than in ages. Table 2 presents the 
ages and [Z/H] values for our sample. 

It is important to note that the age/metallicity anal- 
ysis presented here depends on the assumption of a single 
(luminosity weighted) stellar population, which we feel is 
a reasonable one. The presence of additional populations 
(especially if they are luminous) will bias - perhaps signifi- 
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of the galaxy in the [MgFe] vs H/3 diagram. For instance, 
galaxies placed in the region of 1-3 Gyr have less neighbour- 
ing model points (i.e., lower order) than galaxies older than 
8 Gyr (more neighbouring data points = higher order in- 
terpolation). This is because, as we have stated before, the 
degeneracy of the models is stronger for older systems. Ev- 
ery determination of age and metallicity was double-checked 
against the location in the [MgFe] vs H/3 plane for consis- 
tency of the results. We should always bear in mind that 
non-linear interpolations using distant spaced model points 
as in the case of TMB03 model grids, may result in non- 
stable determinations and there is always the possibility of 
wild oscillations between the tabulated points (see discus- 
sion about interpolation algorithms in Press et al. 1992). 

It is important to notice that the ages and metallicities 
derived here should not be taken as absolute values; there are 
systematic offsets if different index combinations are used, 
although the relative ranking of the galaxies in the sample 
should not be very different (cf. Figure 5) . In the TMB03 sys- 
tem, the average age for the group, field and isolated sample 
(excluding the three Virgo cluster galEixies) is 4.2 Gyr, with 
an average metallicity of -1-0.41 dex. The cluster galaxies in 
our sample also present rather young ages (NGC 4365: 3.6 
Gyr, NGC 4374: 3.8 Gyr, NGC 4754: 5.8 Gyr). Both ellipti- 
cal and SO galaxies in our sample show similar spread in age 
and metallicity in the [MgFe]-H/J diagram with a marginal 
tendency for the bulge of SOs to be slightly younger than 
the elliptical galaxies, as will be shown in Section 4. 

In panels a and b of Figure 6 we compare our determi- 
nations with ages and metallicities in the Terlevich & Forbes 
(2002, TF02) catalogue, for which the authors have compiled 
ages, metallicities and abundances from the highest quality 
data of early-type galaxies in the literature. In TF02, ages 
and [Fe/H] were derived by interpolating Worthey (1994, 
W94) stellar population models in the H/3-[MgFe] plane. To 
carry out this comparison, we have also interpolated ages 
and metallicities using W94 models. 

Figure 6 shows a reasonably good agreement, inside the 
errors, between 15 galaxies in common with TF02 (see also 
Table 1). Our ages are in general slightly younger most prob- 
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Figure 7. The top panel shows the H/3-[MgFe] plane from where Age and [Z/H] were interpolated using the TMB03 model grids. The 
model grids are shown in age steps of 1, 2, 3, 5, 8, 12, 15 Gyr and metallicity steps of -0.33, 0.00, +0.35, +0.67 dex. The derived age 
vs [Z/H] diagram is shown in the bottom panel. Open circles represent elliptical galaxies and filled triangles axe SO-type galaxies in the 
sample. Large circles and the large triangle are isolated galaxies. The three Virgo cluster galaxies in the sample, NGC 4365, NGC 4374, 
NGC 4754, are plotted as square symbols (SO is open square, Es are filled squares). The average la error ellipse is presented at the 
bottom left corner of the age vs [Z/H] diagram. The vectors shown in the error ellipse correspond to 0.078 dex in [Z/H] and 0.17 in log 
Age. 



© 2004 RAS, MNRAS 000, 1-23 



//. Global trends 9 



Table 2. Age, [Z/H] and [a/Fe] from nuclear re/8 extractions. 
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Galaxy names in bold characters are ellipticals. 

The > symbol represents a lower limit estimate: 15.0 Gyr in Age and 0.67 dex in [Z/H]. 
t Derived assuming model tracks of a constant age of 4 Gyr. 



cantly - the resulting age/metallicity estimates. (Bressan et 
al. 1996; Trager et al. 2000a; TOOb). 



4.1 Errors 

It is known that correlated observational errors can mas- 
querade as real trends. Trager et al. (1998) have demon- 
strated the role of index errors in the determination of ages, 



mctallicitics and their correlations with velocity dispersion. 
Furthermore, the Balmer lines used to infer ages and metal- 
licities in e.g. H/3-[MgFe] diagrams, can be severely affected 
by emission and its correction is crucial to the determina- 
tion of galaxy parameters. The general effect of correcting 
H/3 by the presence of emission is to change both the derived 
ages and metallicities. Age is reduced and [Z/H] increased 
([a/Fe] is less affected). It happens so that the informa- 
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tion on H/3 errors is commonly used to separate and select 
the best galaxy data in the literature. The sample used by 
Trager et al. (1998) had a typical error of <5H/3 = 0.191 A. 
A reasonable guide is that H/3 must be accurate to ~ 0.1 A 
to determine reliable ages and mctallicitics. 

Our measured observational errors are not as small as 
Gonzalez (1993, 5H/3 = 0.060 A) or Kuntschner (1998, 6R/3 
= 0.089 A), although these are estimated errors. The average 
H/3-index error in our sample is 0.11 A. Our errors have been 
increased to account for uncertainties in the H/3 emission 
correction: previous to the emission correction our H/3 mean 
error was 0.086 A. 

It is important to emphasize here that in this work we 
focus our attention to conclusions drawn from the relative 
position of galaxies with respect to each other and not on 
absolute values. 

The independent errors in the observables H/3 and 
[MgFo] from Figure 5 translate into correlated errors in the 
age vs metallicity plane. We have calculated the uncertain- 
ties in age and metallicity by interpolating new age and 
[Z/H] values for the points: (-ff/3— error, [MgFe]), {Hl3+error, 
[MgFc]), (H/3, [MgFe]- error) and (H/3, [MgFe]+ error). We 
show in Figure 8 that these four points describe well the 
uncertainties in the central age and [Z/H] determinations. 

We have performed Monte-Carlo simulations of a large 
number of random points inside the la error contour of a 
galaxy in the H/3- [MgFc] plane. The random points were 
then interpolated using the TMB03 model grids; the age 
and [Z/H] values derived are plotted on the right-hand side 
panels of Figure 8. Note that the la error ellipses on the H/3- 
[MgFe] plane were transformed into rather twisted ellipses 
on the age- [Z/H] plane. 

The age and [Z/H] values inter- 
polated for (iJ/3- error, [MgFe]), {H/3+ error, [MgFe]), (H/3, 
[MgFe]— error) and (H/3, [MgFe]+ error) are also shown in 
the right-hand side panels. Although these four points do 
not describe perfectly well the total range of possible age 
and [Z/H] values inside the la error ellipses, they do pro- 
vide us with a sufficient sampling of the uncertainties in- 
volved when transposing errors from the H/3- [MgFe] plane 
to ages or metallicities. We have then used the range of val- 
ues from these four points as the uncertainty in age and 
[Z/H] determinations presented in Table 2. 

We also show in Figure 9 how the error ellipses of Figure 
8 are distributed in the total age- [Z/H] plane. Note that the 
slope of the error ellipses is very close to the direction of 
the spread of the points in the age and metallicity plane 
(we estimate A log age = —1.57 A [Z/H]), indicating that 
correlated errors can masquerade as real trends. 



5 THE [a/FE] RATIO 

According to several authors, the central regions of giant 
E and SO galaxies show an overabundance of a-clomonts 
over Fe, indicated in particular by the steep Mg2 index ra- 
dial gradients (Worthey, Fabcr & Gonzalez 1992; Trager et 
al. 2000a; TOOb). This finding strongly impacts on the the- 
ory of galaxy formation, as super-solar [a/Fc] ratios seem 
to require short star formation time-scales (;$ 1 Gyr) that 
are hardly achieved by current models of hierarchical galaxy 
formation (e.g., Thomas & Kauffmann 1999). 



The TMB03 predictions are identical to TOOb Model 1, 
except that the latter do not include the a-clcmont Ca in the 
a-element group. As pointed out by TMB03, the fractional 
contribution of Ca to the total metallicity is too small (~ 0.1 
per cent) to change the isochrone and SSP characteristics. 

Super-solar [a/Fe] ratios at fixed total metallicity are 
produced by a decrease of the Fe abundance rather than by 
an increase of the a element abundances ([Z/H] = [Fe/H] 
-I- A [a/Fe], where A = - A [Fe/H] /A [a/Fe]). The actual 
atomic abundance ratios of the so-called enhanced elements 
are in fact virtually solar; it is only the Fe-peak elements - 
which contribute ~ 8% to the total amount of metals in the 
Sun - that are depressed (e.g., Trager et al. 2000a, TMB03). 
This result has been interpreted as being due to a fast early 
chemical enrichment, therefore dominated by Type II super- 
novae, in the nuclei of giant E and SO galaxies. Alternative 
explanations, however, include a changing IMF slope and/or 
a changing binary fraction for Type I supernova production. 
Worthey (1998) discusses the possible scenarios thoroughly 
in his paper, and concludes that not all a-elements follow 
the Mg trend. Worthey (1998) even suggests the existence of 
a third chemical enrichment source (perhaps another type of 
SN) other than SN Type I and II, to explain the behaviour 
of N or of the V, Sc, Ti trio in the TiOi,2 bands. Comparing 
[Mg/Fe] ratios from spiral bulges (McWilliam & Rich 1994; 
Jablonka et al. 1996), SOs and elliptical galaxies in the liter- 
ature (Davies et al. 1993; CaroUo & Danziger 1994; Fisher et 
al. 1995, 1996), the data available to date seem to indicate 
that [Mg/Fe] is near zero in most galax;ies of all types with 
velocity dispersion less than ~ 200 km/s (Worthey 1998). 

We will examine below the behaviour of some important 
indices with varying [a/Fe] models. In Figure 10 the indices 
of the a-element group are plotted against the average iron 
index <Fe>. These diagrams are almost degenerate in age 
and metallicity, and so are ideal to identify [a/Fe] variations 
by comparing TMB03 models with [a/Fc] = 0.0 and +0.3. 

Assuming that the difference in central velocity disper- 
sion aa is determined by the depth of the gravitational po- 
tential and ultimately the mass of the galaxy, wo separate 
galaxies according to mass in Figure 10 by assigning differ- 
ent symbols for galaxies with ctq smaller (filled circles) or 
greater (open circles) than 225 km/s. 

The formal [a/Fe] ratios per galax:y were determined 
by interpolating [a/Fe] values from the almost degenerate 
model grids in the Mg2/<Fe> and Mg 6/<Fe> diagrams, 
and averaging the resulting values. The interpolation is per- 
formed between TMB03 models with [a/Fe] = 0.0, -hO.3, 
-1-0.5 and assuming, for simplicity, only model tracks for an 
average age of 4 Gyr (i.e., to make the interpolation eas- 
ier; other model tracks do not alter the results because the 
Mg/Fe plane is strongly degenerate). We present the derived 
[a/Fe] for the galaxies in our sample in Table 2. It is very 
difficult to estimate errors for this determination of [a/Fe], 
but examining the errors in Mg and Fe-indices, we assume 
that on average ± 0.05 dex in [a/Fe] is a very safe upper 
limit estimate of the uncertainty. 

We find that the maximum magnitude of the [Mg/Fe] 
overabundance is around 0.3 dex, for the central re/8 extrac- 
tion. In general, oven galaxies with gq < 225 km/s show Mg 
enhancement with respect to Fe, and the distinction of both 
samples (with ao < 225 km/s and with ao > 225 km/s) is not 
always clear in Figure 10. The NaD and CN indices, trac- 
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Figure 8. Examples of index errors transposed to the age-metallicity diagram. The [MgFe]-H/3 plane is shown on the left-hand side and 
overplotted are the model grids of TMB03. In each plot we show the position of a galaxy and its error bars in the index-index plane. On 
the right-hand side panels we show the results of the model grid interpolation to the age-metallicity plane. The Icr error ellipses were 
obtained from Monte Carlo simulations. 



© 2004 RAS, MNRAS 000, 1-23 



12 G. Denicolo et al. 




Figure 9. Age-metallicity diagram with error ellipses. The error ellipses are the same shown for four galaxies in Figure 8. The slope of 
the error ellipses is nearly the same as the slope of the overall tendency between age and metallicity (we estimate A log age ~ —1.57 A 
[Z/H]), showing that correlated errors can masquerade as real trends. 



ers of Na and mainly N respectively (Tripicco & Bell 1995, 
Worthey 1998, TMB03), show the same type of behaviour 
as Mg with respect to Fe. It appears that Na is more en- 
hanced than Mg, or simply that our measurements are suf- 
fering from NaD absorption by interstellar material. Except 
for a few of the less massive galaxies ( cto < 225 km/s) the 
rest show a-element enhancement (see Figures 10 and 11). 
The Ca4227 and Ca4455 indices barely overlay the models 
for [a/Fe] = 0.0, and their different behaviour with respect 
to the other indices cannot be assigned to [a/Fe] enhance- 
ment. Both indices are not well calibrated for galaxies, and 
not convincingly either for globular clusters (Maraston et al. 
2003). Ca4227 is sensitive mainly to N, and its Ca-element 
dependence is yet to be shown. As for Ca4455 - which is 
insensitive to Ca abundance (Tripicco & Bell 1995), Fe and 
Cr seem to be the main contributors. The TiO model indices 
are almost degenerate in the [a/Fe] sensitivity, and are badly 
calibrated even for globular clusters in Puzia et al. (2002). 

In summary, if large [a/Fe] ratios require short star for- 



mation time scales for the galaxies, we can say from Figure 
10 that, in our sample of mostly low-density environment 
galaxies, massive early- type galaxies (high ctq) do not ob- 
viously have more rapid star formation (high [a/Fe]) than 
small (i.e., less massive) early-types, as predicted by hierar- 
chical galaxy formation models. 

Figure 11 shows the relations of the [a/Fe] ratio 
with age, metallicity, velocity dispersion and magnitude 
per galaxy. There is a trend in that older galaxies are 
strongly overabundant whereas younger galaxies approach 
solar abundance ratios (Figure 11-a). A Spearman rank cor- 
relation test results in p = 0.38 with more than 90% prob- 
ability that a correlation is present. It is necessary to note, 
nevertheless, that part of this correlation may be due to the 
correlated errors discussed before. 

The median [a/Fe] ratio for the 50 Es is 0.19 ± 0.08 dex, 
and for the 30 SOs it is 0.11 ± 0.11 dex, the uncertainties 
are the rms values. 
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Figure 10. Index tracers of a-elements vs <Fe>. Overplotted are the a-enhanced SSPs by TMB03 for [a:/Fe]=0.0 (solid lines) and 
[a/Fe]=+0.3 (dashed lines). The symbols separate galaxies with ao smaller or greater than 225 km/s, as indicated in the top-left panel. 
The outlier galaxy on the bottom-left of all panels is NGC 3156. 



Trager et al. (2000b, TOOb), analysing the sample of 
naainly "field" elliptical galaxies from Gonzalez (1993), 
claimed that [a/Fe] depends only on uo and not on age. 

However, Figure 11-a suggests a SNe II/SNc la enhance- 
ment ratio higher in the older galaxies of our sample; in 



other words, there is a hint that older galaxies had short 
star-formation time scales. On the other hand, there is no 
clear correlation of [a/Fe] with metallicity in Figure 11-b. 

TOOb found that [a/Fe] and cto are highly correlated 
(linear) quantities. Their proposed linear fit to the relation 
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Figure 11. Correlations of the enhanced element ratio [a/Fe] witli (a) age, (b) metallicity, (c) velocity dispersion and (e) absolute 
magnitude. Panel (d) shows the scatter about the solid line fit from panel (c) and the dotted lines represent the Icr scatter (A[a/Fe]=0.079 
dex. The median [a/Fe] ratio of the sample is 0.17 dex. Circles reijrcsent elliptical galaxies and triangles are SO-type galaxies in the 
sample. Large circles and triangles are isolated galaxies, medium size denote field, and the smallest symbols represent group galaxies. 
The three Virgo cluster galaxies in the sample, NGC 4365, NGC 4374, NGC 4754, are plotted as square symbols (open square denotes 
SO, filled square denotes E). The error bars represent the mean uncertainties in log Age (0.17 dex), [Z/H] (0.078 dex), log ao (0.024 dex) 
and [a/Fe], for which we assume the error estimate of 0.05 dex. 



is shown by the dashed hno in Figure 11-c. Examining panel 
(c) one cannot safely conclude that the previously discussed 
trend in [a/Fe] is as prominent: we find the linear correlation 
coefficient R = 0.63. Most of the spread in Figure 11-d seems 
consistent with the errors, thus, for a given ao most (but 
not all) galaxies have a similar [a/Fe] within the errors. The 



existence of [a/Fe] ratio variations at a given central velocity 
dispersion would imply that the star-formation history of 
old galaxies at a given cro is not entirely homogeneous. The 
[a/Fc]-a() fit to our sample is shown by the solid line in 
Figure 11-c (slope: 1.7). We have used the bisector (solid 
line) of the regressions of [a/Fe] on ao and ao on [a/Fe] 
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Table 3. Principal component analysis: the hyperplane. 
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Figure 12. The edge-on view of the mctallicity-agc-cr plane. The 
dashed line in the top panel is the fit to our data points and 
defines the plane as [Z/H] = 0.99 log ctq - 0.46 log Age - 1.60. 
The dashed line in the bottom one is the fit in Trager et al. 
(2000b). Vectors of 0.17 in log Age, 0.078 dex in [Z/H] and 0.024 
in log (TO represent the average of the overall uncertainties in our 
sample. The symbols are the same as in Figure 11. 



(represented by the two dotted lines in panel c). Panel (d) 
shows the residuals of the [a/Fej-ao relation (with respect 
to the solid line of panel c); the la scatter of the relation 
is A [a/Fe] = 0.079 dex, and is represented by the dotted 
lines in panel (d). 

In Figure 11-c we present the weak relation between 
[a/Fe] and absolute blue magnitude. There is quite a large 
range of magnitudes for oversolar [a/Fe] values. 

In sum these results reinforce the picture of very old 
and massive galaxies having formed in relatively short time 
scales with small scatter in their [a/Fe] values. 



6 [Z/H], AGE, (70 AND [a/FE] HYPERPLANE 

In a very interesting paper, Bressan et al. (1996) explored 
the four dimensional space which involved H/3, [MgFe], 
(1550-V) and log ctq in elliptical galaxies. Forbes & Ponman 
(1999) were among the first to explore the relation between 
(Jo and stellar population ago derived from model compari- 
son. Later, TOOb performed a principal component analysis 
(PCA; other references to PCA: Whitney 1983, Ronen et al. 
1999) on four variables, namely log ao, log Age, [Fe/H] and 
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[E/Fe]; i.e. velocity dispersion, ago, metallicity and element- 
enhanced ratio which they called E/Fe (here wc adopt the 
notation a/Fe). TOOb found that most of the variance (91 
percent) in their sample of local ellipticals can be explained 
with just two parameters, i.e., the sample is basically bi- 
paranictric in the four-dimensional space ao, Age, [Fe/H] 
and [E/Fc], in other words, the galaxies are located in a 
hyperplane in this four-dimensional space that they named, 
the metallicity hyperplane. 

Given the importance of this finding, wo have performed 
a PCA in our sample data. The results arc presented in Table 
3 and Figure 12. 

It is interesting to note from an inspection of the co- 
variance matrix that while log a is significantly correlated 
to the other three variables (log Age, [Z/H] and [a/Fe]), 
both log Age and [a/Fe] are not correlated to [Z/H] in our 
sample. A similar conclusion regarding the behaviour of log 
(To, [Z/H] and log Age can be obtained from Figure 13 where 
the panels show the sample segregated according to ao. 

The least square fit to our hyperplane is: 

[Z/H] = 0.99(±0.15)Zoff<To-0.46(±0.08)/offApe-1.60(±0.40)(l) 

In hnear terms, this equation can be translated to the 
approximate law: 

-0.5 



Z oc ao ■ Age 



(2) 



If we assume constant values for ao and plot eq. 1 in 

the Age-metallicity plane, we find the interesting relations 
shown in Figure 13. The existence of the Z-plane indicates 
that there exists an age-metallicity relation for each value 
of ao- This relation implies, as shown in Figure 13, that at 
a fixed (To, metal- rich galaxies have a younger stellar pop- 
ulation in their central regions. TOOb have found the pro- 
jection of the hyperplane to be [Z/H] = 0.76(±0.13) logiro 
- 0.73(±0.06) log Age - 0.87(±0.30), where Age is always 
taken in Gyr. 

From the correlation matrix, the strongest effect seems 
to be a change in log age (and not in [Z/H]) in the sense 
that age increases going from small ao ellipticals to large 
ao ones. Even though there are several lower limits in the 
right panel of Figure 13 (massive galaxies that fall outside 
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Figure 13. Age-metallicity plane for constant ctq- Circles represent elliptical galaxies and triangles are SO-type galaxies in the sample. 
Large circles and triangles are isolated galaxies, medium size denote field, and the smallest symbols represent group galaxies. The three 
Virgo cluster galaxies in the sample, NGC 4365, NGC 4374, NGC 4754, are plotted as square symbols. The panels separate the galaxies 
in bins of <To; the dashed lines are fits from eq. 1, assuming constant ao values as labelled. 



the model grids) the effect is visible between the left and 
central panels. 

This is an interesting result, and cannot be a reflection 
of the well known Mg-a relation. If that was the case one 
would expect a stronger correlation between ao and metal- 
licity, given that Mg is the dominant index in the [Z/H] 
estimate in this work using TMB03 models. 

Nevertheless, one must be aware in dealing with projec- 
tions of hyperspaces that the apparent correlation may only 
bo the product of the distribution of points in the hyper- 
space projected onto a particular plane. 

To answer this point we performed a Principal Compo- 
nents Analysis, the results of which are given in Table 3. The 

first two principal components contain 84 percent of the to- 
tal variance indicating that the ellipticals in our sample are 
distributed in a flattened volume in the 4-D space. This is 
further illustrated in Figure 14 where we show the face-on 
and edge-on views of the projected variables. To simplify 
comparisons with TOOb we have used similar orientation in 
our diagrams. It is interesting how similar are our results 
and those of TOOb, strengthening the reality of the metal- 
licity hyperplane. 



Summarising, metallicity is a linear function of both 
log Age and log ao, i.e., Z oc ao-Age~°'^ , which defines 
the so called "Z-plane". Following the fit we derived for 
our sample, the contours of constant ao have the slope 
Alog^5e«-2.2A[Z/H]. This is close to the "3/2 relation" 
of Worthey (1994), which expresses trajectories in the log 
Age - [Z/H] space along which colours and line-strengths 
remain roughly constant. Thus, the Z-plane correlation pre- 
dicts that line-strengths should be constant along trajecto- 
ries of constant tJo in the Z-plane. 

The enhancement ratio, [a/Fe], is a linear function of 
log (To, increasing towards high-cro galaxies, and should be 
also responsible for some of the scatter in Figure 12. 



7 DISCUSSION 

Kuntschner (1998) found that the Fornax cluster elliptical 
galaxies appear to be roughly coeval with some evidence for 
slightly younger ages for the most metal-rich ones. Using 
Lick/IDS indices and Worthey (1994) SSP models, he esti- 
mated the ages to be around 8 Gyr, and the galaxies showed 
a sequence of metallicities from -0.1 to -1-0.35 in [Fe/H]. This 
result is consistent with the conventional view of old, coeval 
elliptical galaxies where the metallicity scales roughly with 
the size of the galaxy. However, the lenticular galaxies in 
Fornax have luminosity weighted ages that arc smaller than 
those of the ellipticals, spanning from less than 2 Gyr to 

8 Gyr, and also covering a large range of metallicities from 
-0.5 < [Fe/H] < -1-0.5. For the purpose of comparison with 
fiterature data, Kuntschner (1998) also observed the galaxy 
NGC 3379. The age of NGC 3379 in Kuntschner's analysis 
is ~ 8.4 Gyr, whereas we have estimated ~ 10.0 Gyr (with 
TMB03 models) . We feel that comparing age estimates from 
different authors (which may have used different models) is 
not a reliable procedure; to check for systematic differences 
between cluster and field environment we prefer to compare 
the distribution of the galaxies in the [MgFe]-H/J plane with 
the Fornax sample distribution, as both are in the Lick/IDS 
index system. 

In Figure 15 we show a comparison between H/3 and 
[MgFe] indices for the Fornax and OAGH samples. The bot- 
tom panel of Figure 15 shows only SO galaxies in both sam- 
ples, and the right-hand panel shows only elliptical galaxies. 

We estimate that the average age for our sample of 
group, field and isolated galaxies (5.8 Gyr) is younger than 
Fornax cluster ellipticals (~ 9-10 Gyr). Secondly, both Es 
and SOs in our sample cover a large range of ages and metal- 
licities. Although our galaxies span wider ranges in age and 
[Z/H] than the Fornax galaxies, the Es and SOs in the OAGH 
sample appear to have on average more similar stellar pop- 
ulations than the very distinct behaviour of Es and SOs in 
Fornax. 

In Section 4 we have estimated that SOs present a 
younger mean age than Es. It is important to emphasize 
here that the line-strength indices reflect only the integrated. 
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Figure 14. Principal Component Analysis results on the OAGH sample. The position of the galaxies in the diagrams is denoted by 
their names. The top panels represent the face-on (right) and edge-on (left) projections of the log Age-[Z/}l]-log <7-[a/Fe] hyperplane. 
Projections of the four basic variables are shown as arrows in the direction of increase (for log Age, this arrow points in the direction 
of older galaxies). Velocity dispersion and enhancement ratio dominate the first principal component (PCI), while age and metallicity 
dominate the second (PC2). The third (PCS) and fourth (PC4) principal components contribute less than 12% to the overall variance 
in the logAge-[Z/H]-log a-[a/Fe] space. 



luminosity weighted properties in a galaxy. As young pop- 
ulations are more luminous per unit mass than old ones, a 
small (in mass) young population can dramatically change 
the strength of indices, in particular H/3. de Jong & Davies 
(1997) and Trager (1997) have demonstrated this and also 
showed that the strong H/3 galaxies in Gonzalez (1993) sam- 
ple tend to have disky isophotes. de Jong & Davies (1997) 
suggested that ongoing star-formation might be associated 
with the presence of a disk. In the OAGH sample, ellipti- 
cal galaxies that present young ages and hints of a stellar 
disk are NGC 315 (Pcleticr et al. 1990; our age: 2.5 Gyr), 
NGC 584 (Michard & Marchal 1994; 3.8 Gyr), NGC 720 
(Goudfrooij et al. 1994; 3.0 Gyr) , NGC 821 (Ravindranath 



et al. 2001; 4.0 Gyr), NGC 3377 (Scorza 1993; 3.6 Gyr), 
NGC 3610 (Scorza & Bender 1990; 1.8 Gyr), NGC 3613 
(Goudfrooij et al. 1994; 5.6 Gyr), NGC 3640 (Goudfrooij 
et al. 1994; 2.5 Gyr), NGC 4365 (Forbes et al. 1995; 3.6 
Gyr), NGC 4374 (Zirbel & Baum 1998; 3.8 Gyr), NGC 5322 
(Goudfrooij et al. 1994; 2.4 Gyr), NGC 5845 (Reid et al. 
1994; 2.5 Gyr), i.e., 60% of elliptical galaxies with ages 
younger than 4 Gyr show photometric evidence of a stel- 
lar disk component in their morphologies. 

NGC 3156 has a very distinct behaviour with respect 
to the bulk of the sample. As we have noted in some plots, 
this is the youngest and most metal-poor elliptical galaxy 
in the sample (NGC 5864 shows equivalently low metallicity 
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Figure 15. Comparison between the distribution in the H/3-[MgFc] jjUinc of Fornax and OAGH samples. The model grids are from 
TMB03, with [a/Fe] = 0.0, in age steps of 1, 2, .3, 5, 8, 12, 15 Gyr and nictallicity steps of -0.33, 0.00, +0.35 and +0.67 dex. The bottom 
panel shows only SO galaxies in both samples, and the right-hand panel shows only elliptical galaxies. 



but with larger error bars, and substantially older age). Ac- 
cording to the calculations of Poggianti & Barbaro (1997) 
strong Balmer absorption is created mainly at the main se- 
quence turn off and cannot originate from other evolutionary 
stages (at least not for relatively young systems) . Note that 
we are referring to strong Balmcr absorption in a galaxy 
with a young stellar population. In many cases. Horizontal 
Branch (HB) stars can affect Balmer lines rather severely 
(see Maraston et al. 2003 for TMB03 model calibration of 
the Balmer lines). However, in the case of NGC 3156 it would 
be necessary to have a significant population of HB stars to 
produce such a young age. In fact, Maraston et al. (2003) 
have determined rather old ages (~ 8 Gyr) for two globu- 
lar clusters NGC 6441 and NGC 6388 which are known to 
exhibit an important tail of warm HB stars. Furthermore, 
Burstein et al. (1988) have shown that the ultraviolet excess 
in early-type galaxies (which is now known to be associated 
with HB stars; O'Connell et al. 1999) correlates well with 
Mg2 (NGC 3156 shows weak Mg2 absorption [cf. Paper I], 
and poor mctallicity) . Hence wc believe that the young age 
estimated from the model predictions for NGC 3156 is reli- 
able. This is a disk dominated SO galaxy, with a dusty and 
faint bulge, located in a group and shows no obvious inter- 
action features. It is a galaxy that shows a high V/ao value 
(=1.2; Paper I). 

The relations between Age and metallicity [Z/H] against 



velocity dispersion and absolute blue magnitude (calculated 
from apparent magnitudes and our determinations of radial 
velocity given in Paper I, and assuming Ho = 75 km/s/Mpc) 
are shown in Figure 16. We observe a trend in that the old- 
est and more metal rich galaxies tend to have higher velocity 
dispersions. Caldwell, Rose & Concannon (2003) and Trager 
et al. (2000a,b) claim that there is a trend for lower ao (or 
less massive elliptical) galaxies to have younger ages, con- 
trary to what is found by Kuntschner (2000) and TF02 with 
their samples. Caldwell et al. (2003) also note a modest (if 
any) trend in [Z/H] with a. Ours and these previous results, 
however, are not conclusive enough to place key constraints 
on age distributions predicted by numerical simulations of 
galaxy formation, but we refer the reader to a discussion of 
this point in TF02. 

In Section 5 we have seen that there is a marginal ten- 
dency for massive galaxies to show larger [a/Fe] values than 
less massive ones. Worthcy (1992,1998) detail several pro- 
cesses that can produce an Q-elemcnt enhancement as a 
function of galaxy mass. They are mostly related to the 
different contributions towards the chemical enrichment of 
galaxies by supernovac types la and II. SN la (with progen- 
itors less massive than 5-10 M0) release Fe-pcak elements 
after a time delay of about a few times 10* years. Type II 
supernovae, on the other hand, originate from more massive 
stars and contribute earlier to the abundance of a-elements 
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Figure 16. Panels a and b present the central velocity disper- 
sion against Age and metallicity determinations for the re/8 aper- 
ture extraction. Panels c and d show the absolute blue magnitude 
against Age and [Z/H]. The symbols are the same as in Fig. 11. 



and Fe. The ratio of SNe type II to type la, which can vary 
from massive to less massive galaxies, somehow determines 
the [q/Fc] abundances. There are several scenarios proposed 
for the variation of type Il/type la ratio (e.g. Arimoto & 
Yoshii 1987) including winds, amount of binaries, even IMF 
variations due to environmental causes. Different time-scales 
for the formation of SNe Types II and la, allow the predic- 
tion of a correlation of Mg strength with dynamical crossing 
time (Worthey 1998) or, if including the effects of winds 
(Arimoto & Yoshii 1987) with the potential well of the 
galaxy. An environmentally controlled IMF predicts good 
correlation with a (Paper I). While Mg2 versus cr is a well 
known tight correlation, <Fc> versus a is weaker. This sug- 
gests that the basic relation is between the Mg abundance 
and the velocity dispersion. Contrary to Worthey (1998), 
who finds that [Mg/Fe] vs. a does not depend on galaxy 
type, supporting the statement above, we find that ellipti- 
cals have [a/Fe] ratios slightly higher than SOs, on average. 

Other indices that correlate well with a are NaD and 
CN (Worthey 1998). Matteucci & Padovani (1993) claim 
that a change in the IMF slope of 0.4 could produce a 0.3 dex 
change in a-element abundances. It remains to be proved 
that such a change in the IMF slope could be consistent 
with observations. 

The metallicity of the interstellar medium can provide 
a constraint to decide between the time formation difference 
hypothesis and the IMF variation one (Worthey 1998). We 
have to resort to X-ray abundance determinations where we 
encounter difficulties with modelling the Fe atom (Loewen- 
stein & Mushotzky 1996; Arimoto et al. 1997) although 
progress has been achieved in that direction (Matsushita 
et al. 2000). 

Finally, in the previous Sections we have seen that, in 
general, high and low density environment galaxies (the for- 



mer in groups and the latter field or isolated) have sim- 
ilar distribution in the Age- [Z/H], H/3-[MgFe], H7-[MgFe] 
planes, thus suggesting that dynamical interactions are not a 
necessary driver for the observed variations. If this is true, it 
confirms Longhctti et al. (1998a,b, 1999, 2000) who reached 
the same conclusion by examining the line-strength indices 
of 51 early- type galaxies (21 shell galaxies and 30 members 
of pairs displaying clear signatures of interaction) and com- 
paring them with the reference sample of "normal" galaxies 
by Gonzalez (1993). 

One should be prudent not to overinterpret these re- 
sults, as the truly isolated galaxy statistics in our sample is 
limited by small number effects. We coxild in fact claim to 
be seeing a different distribution on the H/3,[MgFe] plane for 
Fornax as compared with our isolated/field/group galaxies 
in Figure 15. 

Our new data seem to be showing a trend for massive 
(isolated/field/group) galaxies to be older (Figure 13 and 
16-a) and display greater a/Fe values (Figure 11-c). 



8 A BROAD PICTURE OF GALACTIC 
CHEMICAL EVOLUTION 

Understanding the broad picture of galactic chemical evolu- 
tion will require us to firm up the links between the chemical 
abundances of predominantly active star-forming systems 
(i.e. spirals) and predominantly quiescent ones (i.e. ellipti- 
cals, spheroids, lenticulars and bulges). 

While a common elemental yardstick may not exist be- 
cause of the different elements which we observe directly in 
each galaxy type, it may be possible to tie the two types to- 
gether abundances- wise by observing elements in each which 
share the same nucleosynthesis production site. An example 
might be oxygen and magnesium. In external spirals oxygen 
is taken as the metallicity gauge primarily because of its 
observability. Magnesium, which, like oxygen, is primarily 
produced in massive stars (Nomoto et al. 1997a,b) may be 
measurable directly through a calibrated Mg2 index. Then 
oxygen and magnesium might be linked by assuming a "cos- 
mic" Mg/O ratio calibrated locally. Some features in the 
optical spectra of elliptical galEixies that contain oxygen sig- 
natures are G4300 and C24668, in the near-infrared there is 
also the 2.3 /im CO feature. As synthetic spectra and stellar 
abundances grow more precise, the contribution of specific 
elements to the indices might be measurable, leading to a 
much more clear understanding of chemical enrichment and 
galaxy formation. 

In what follows we have used the Hq line information 
in our spectra of early-type galaxies to estimate the oxygen 
abundance using the N2 ratio as in Denicolo, Terlevich & 
Terlevich (2002). We have derived the metallicity [Z/H] for 
our ellipticals and SOs in Section 4 (this metallicity comes 
mainly from Mg and Fe abundances), so it is an interesting 
exercise to compare both [Z/H] metallicity and O abundance 
estimate for early-type galaxies. 

The N2 calibrator is defined as the logarithmic ratio of 
the emission lines [Nil] A6584 to Ha. In the elliptical galax- 
ies, we have measured [N ll] and Ha as indices. The N2 is the 
logarithmic ratio of fiuxes and not equivalent widths, so we 
have computed fluxes from the index definitions as stated 
in eq. 11 of Paper I. We have only used [N ii] and Ha in- 
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Figure 17. Oxygen abundance vs metallicity [Z/H]. The oxygen abundance was derived from the [Nil] and Ha index measurements in 

the spectra of our early- type galaxies using the N2 calibrator. The [Z/H] metallicity is the result of the interpolation of Thomas et al. 
(2003) SSP models for measurements in galaxy spectra within the re/8 aperture. Solid circles indicate elliptical galaxies, open circles are 
SOs. The size of the points is associated with the galaxy age (ages are from Table 2). The dotted line indicates when 12+log(0/H) is 
equal to [Z/H]. The dashed line represents 12-|-log(0/H) when N2 is null. 



dices with negative values (i.e. lines in emission). The errors 
in the [N ll] and Ha indices arc derived from repeated mea- 
surements. The 12-|-log(0/H) abundances were derived from 
12+log(0/H) = 9.23(±0.02) + 0.79(±0.03)xN2 (Denicolo, 
Terlevich & Terlevich 2002). In Table 4 we show the oxygen 
abundance and [Z/H] metallicity, for the galaxies with [Nil] 
and Ha indices in emission. 

In Figure 17 we present the 12-Mog(0/H) vs [Z/H] 
diagram. The plot supports substantial oversolar abun- 
dances for the interstellar medium (ISM) as measured from 
12-Mog(0/H). If we take the Solar value of 12+log(O/H)0 = 
8.87 (Grevcsse et al. 1996), and by definition [Z/H]0 = 0.0, 
we draw the dotted line in Figure 17 where 12-|-log(0/H) 
= [Z/H]. The distribution of points indicates that most 
galaxies have Z(ISM) > Z (stars), where we use here "Z" 
to loosely describe metallicity; Z(stars) is taken from [Z/H] 
and Z(ISM) is inferred from 12-hlog(0/H). Having Z(ISM) 
> Z(stars) is consistent with our understanding of star for- 
mation and evolution: stars cannot have abundances that 
arc larger than those of the gas from which they formed. 
The plot in Figure 17 suggests that O/H abundance in the 
nuclear ISM is similar for all of our sample galaxies, irrespec- 
tive of the luminosity averaged ages and stellar abundances 



inside re/8. We regard this result as significant. Denicolo, 
Terlevich & Terlevich (2002) showed that the relation N2 
to O/H exists for either black-body or power-law photoion- 
ization. There axe several sources of scatter in that relation: 
variations of the photoionization paxameter, shape of the 
ionizing continuum, differences in secondary production of 
N, to site a few. What is important from the present result 
is the apparent lack of scatter suggesting a uniformity of 
the central ISM composition irrespective of the mass of the 
galaxy. If this result is confirmed with a larger sample of 
gala^xies, one of the consequences of it would be to rule out 
galaxy mass related IMF variations. 



In other words, we can interpret it as if chemical evo- 
lution has finished in these galaxies, i.e., most stars have 
already ejected metals to the ISM and, in general, any cur- 
rently ongoing star formation in these galaxies should occur 
on a small scale. This is an interesting idea, and could be 
what distinguishes most of the chemical characteristics from 
early-type to late-type galaxies. 



© 2004 RAS, MNRAS 000, 1-23 



//. Global trends 21 



Table 4. Oxygen abundance and [Z/H] for early-type galaxies. 



Galaxy 12+log(0/H) error 


[Z/H] 


EW{Ha) 


er 


(Ha) 


EW([Nn] ) 


cr([Nn] ) 


NGC 584 


9, 


.376 


0, 


,146 


0.435 


-0, 


,244 


0, 


,081 


-0, 


,287 


0, 


,018 


NGC 821 


9, 


.318 


0, 


,196 


0.480 


-0, 


,316 


0, 


,093 


-0, 


,226 


0, 


,077 


NGC 1407 


9, 


.242 


0, 


,075 


0.690 


-0, 


,531 


0, 


,085 


-0, 


,125 


0, 


,008 


NGC 1600 


9, 


.263 


0, 


,198 


0.690 


-0, 


,545 


0, 


,091 


-0, 


,233 


0, 


,099 


NGC 1700 


9, 


.376 


0, 


,091 


0.498 


-0, 


,389 


0, 


,081 


-0, 


,518 


0, 


,013 


NGC 1726 


9, 


.277 


0. 


,118 


0.362 


-0, 


,800 


0, 


,166 


-0, 


,868 


0, 


.153 


NGC 2300 


9 


.498 


0, 


,295 


0.530 


-0, 


,155 


0, 


,093 


-0, 


,380 


0, 


.122 


NGC 2418 


9 


.164 


0, 


,037 


0.569 


-1, 


,118 


0, 


,091 


-1, 


,356 


0, 


.037 


NGC 2768 


9, 


.032 


0, 


,036 


0.381 


-1, 


,317 


0, 


,097 


-1, 


,594 


0, 


,065 


NGC 2872 


9, 


.312 


0, 


,289 


0.350 


-0, 


,256 


0, 


,115 


-0, 


,162 


0, 


,080 


NGC 3156 


9, 


.321 


0, 


,095 


-0.170 


-0, 


,702 


0, 


,153 


-0, 


,964 


0, 


,022 


NGC 3245 


8, 


.168 


0, 


,038 


0.318 


-1, 


,806 


0, 


,083 


-1, 


,643 


0, 


,121 


NGC 3379 


9, 


.276 


0, 


,096 


0.150 


-0, 


,435 


0, 


,085 


-0, 


,198 


0, 


,020 


NGC 3599 


9, 


.400 


0, 


,078 


0.422 


-0, 


,665 


0, 


,107 


-1, 


,602 


0, 


,129 


NGC 3607 


9, 


.217 


0. 


,034 


0.690 


-1, 


,022 


0, 


,080 


-2, 


,206 


0, 


.020 


NGC 3608 


9, 


.284 


0, 


,136 


0.357 


-0, 


,389 


0, 


,087 


-0, 


,190 


0, 


.042 


NGC 3665 


8, 


.401 


0, 


,024 


0.690 


-2, 


,122 


0, 


,083 


-1, 


,150 


0, 


.046 


NGC 3941 


9, 


.289 


0. 


,063 


0.585 


-0, 


,663 


0, 


,084 


-0, 


,560 


0, 


.040 


NGC 4125 


8, 


.956 


0, 


,034 


0.317 


-1, 


,264 


0, 


,087 


-2, 


,644 


0, 


,098 


NGC 4261 


9 


.480 


0, 


,105 


0.437 


-0, 


,493 


0, 


,115 


-1, 


,438 





,102 


NGC 4374 


9, 


.034 


0, 


,037 


0.505 


-1, 


,308 


0, 


,107 


-1, 


,633 


0, 


,040 


NGC 4550 


9, 


.196 


0, 


,135 


0.290 


-1, 


,048 


0, 


,136 


-0, 


,268 


0, 


,076 


NGC 5322 


9, 


.372 


0, 


,092 


0.690 


-0, 


,391 


0, 


,082 


-0, 


,503 


0, 


,015 


NGC 5353 


9, 


.207 


0, 


,049 


0.690 


-1, 


,029 


0, 


,090 


-1, 


,417 


0, 


.102 


NGC 5354 


9, 


.256 


0, 


,266 


0.041 


-0, 


,464 


0, 


,086 


-0, 


,144 


0, 


.084 


NGC 5444 


9, 


.301 


0. 


,132 


0.192 


-0, 


,539 


0, 


,094 


-0, 


,418 


0, 


.104 


NGC 5813 


9, 


.275 


0. 


,095 


0.271 


-0, 


,784 


0, 


,092 


-0, 


,773 


0, 


.143 


NGC 5845 


9, 


.310 


0, 


,118 


0.690 


-0, 


,418 


0, 


,095 


-0, 


,305 


0, 


,045 


NGC 5846 


9, 


.181 


0, 


,039 


0.337 


-1, 


,075 


0, 


,086 


-1, 


,824 


0, 


,078 


NGC 6172 


9, 


.377 


0, 


,058 


0.250 


-0, 


,729 


0, 


,092 


-1, 


,742 


0, 


,082 


NGC 7302 


9, 


.610 


0, 


,695 


0.332 


-0, 


,070 


0, 


,111 


-0, 


,286 


0, 


,060 


NGC 7585 


9, 


.523 


0, 


,102 


0.690 


-0, 


,371 


0, 


,082 


-1, 


,186 


0, 


,091 


NGC 7619 


9, 


.264 


0, 


,090 


0.690 


-0, 


,635 


0, 


,084 


-0, 


,334 


0, 


.053 



9 CONCLUSIONS 

We have applied the stellar population model of Thomas, 
Maraston & Bender (2003) to estimate ages, metallicities 
and abundance ratios in a sample of 83 early-type galaxies 
mostly located in low density regions. 

The stellar population properties derived for each 
galaxy correspond to the nuclear re/8 aperture extraction. 
Both ellipticals and SOs in the sample show a large spread 
in ages (from around 1 Gyr to 15 Gyr, but many galaxies 
also fall out of the model grids) and metallicities (mostly 
from 0.0 to 0.67 dex in [Z/H]). On average, SO galaxies are 
slightly younger than the elliptical galaxies by almost 3 Gyr, 
and present higher [Z/H] than ellipticals by approximately 
0.10 dex. The average age found for the ellipticals is 5.8 ± 
0.6 Gyr and the average metallicity is +0.37 ± 0.03 dex. For 
SOs, the average age is 3.0 ± 0.6 Gyr and <[Z/H]> = 0.53 ± 
0.04 dex. These averages were performed excluding the three 
Virgo cluster galaxies in the sample, and the uncertainties 
shown are errors of the mean. 

It seems that our elliptical galaxies are in general 
3-4 Gyr younger than E galaxies in the Fornax cluster 
(Kuntschner 1998, 2000; see Figure 15). In the Fornax clus- 
ter the E and SO galaxies are easily distinguished in terms of 
age and metallicity, whereas the Es and SOs from low-density 



regions appear to have more homogeneous stellar population 

properties; i.e. the differences in age and metallicity are not 
so noticeable for Es and SOs of low-density regions than for 
Fornax. 

We find that galaxies lie in a plane defined by 
[Z/H]=0.991ogao-0.461ogAsie-1.60, or in linear terms 
Z(xao-Age~'^'^ . In this work we are adding more young galax- 
ies to the coverage of the hyperplane defined in Trager et al. 
(2000b). 

Analysing [a/Fc] ratios, wc can cast some light on 
the star formation timescale of the galaxies. More mas- 
sive (larger <to) and older galaxies present, on average, large 
[a/Fe] values, and therefore, must have undergone shorter 
star-formation timescales. 

Comparing group against field/isolated galaxies, it is 
still not clear that environment plays an important role in 
determining their stellar population history. In particular, 
our isolated galaxies show ages differing by more than 8 
Gyr. However, we note that the cluster elliptical galaxies of 
Kuntschner (2000) are more uniform in H/3 (and age) than 
the group, field and isolated galaxies of our sample. 

Finally, we attempted to draw a broad picture of galac- 
tic chemical evolution by comparing [Z/H] metallicities 
derived from SSP models (i.e. stellar metallicity) against 
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log(0/H) abundances obtained from the N2 ratio of [Nil] 
A6584 to Hq emission lines (i.e. ISM metallicity). The dis- 
tribution of points indicates that most galaxies have Z(ISM) 
> Z(stars), which is consistent with chemical evolution ex- 
pectations. Wo interpret it as if chemical evolution has fin- 
ished in early- type galaxies, i.e., most stars have already 
ejected metals to the ISM and, in general, any currently on- 
going star formation episode in these galaxies should be on 
a small scale. 
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